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The aim of this study is to extend the scope and applicability of the level-crossing method to 
discrete-time stochastic processes and generalize it to enable us to study multiple discrete-time 
stochastic processes. In previous versions of the level-crossing method, problems with it correspond 
to the fact that this method had been developed for analyzing a continuous-time process or at most a 
multiple continuous-time process in an individual manner. However, since all empirical processes are 
discrete in time, the already-established level-crossing method may not prove adequate for studying 
empirical processes. Beyond this, due to the fact that most empirical processes are coupled; their 
individual study could lead to vague results. To achieve the objectives of this study, we first find an 
analytical expression for the average frequency of crossing a level in a discrete-time process, giving 
the measure of the time experienced for two consecutive crossings named as the “waiting time”. We 
then introduce the generalized level-crossing method by which the consideration of coupling between 
the components of a multiple process becomes possible. Finally, we provide an analytic solution 
when the components of a multiple stochastic process are independent Gaussian white noises. The 
comparison of the results obtained for coupled and uncoupled processes measures the strength and 
efficiency of the coupling, justifying our model and analysis. The advantage of the proposed method 
is its sensitivity to the slightest coupling and shortest correlation length. 

PACS numbers: 05.45.Tp,02.50.-r,05.10.-a 


I. INTRODUCTION 

A. Motivation 

The level-crossing method has a long history in 
studying stationary stochastic processes 0 , 0 - The 
matter of interest in the level crossing method is to 
extract the statistical properties of a stationary process 
by its upcrossings through a specific level. This is very 
important from a mathematical point of view, since it 
is closely related to the problem of extremes in random 
processes 0. In the context of the present study, a more 
hands-on practical aspect of level crossing named the 
“waiting time” is of interest. In the level-crossing method, 
the average time between two consecutive crossings of 
the level is a measure of the “waiting time”. The concept 
of waiting time has been studied in various fields, e.g., 
statistical physics 0,0 , material physics quantum 

physics [§, 10, solar physics [Til . fluid dynamics 
|H [l4| , and finance 0 {lTj . 

It should be noted that prior to the present study the 
level-crossing method has been mostly implemented for 
studying continuous-time processes, paying less attention 
on discrete-time processes either for individual or coupled 
processes. Now two questions may arise. Firstly, why 
does studying level crossing for discrete-time processes 
matter? secondly, why are we required to study coupled 


processes? The significance of studying level crossing 
for discrete-time processes comes from the fact that all 
empirical processes are discrete in time. This importance 
has been amplified by the inconsistency observed 
between the level-crossing results for continuous- and 
discrete-time processes. Therefore, we need to study 
the level-crossing method for discrete-time processes 
before applying it to empirical data. To address the 
second question, we should note that the level-crossing 
method had been introduced for analyzing a single 
process. In other words, it works fine for the individual 
study of processes. However, we sometimes need to 
study multiple stochastic processes, where some sort of 
coupling or coexistence between processes exist. This 
motivated us to generalize the level-crossing method 
for studying multiple discrete-time processes. The 
generalized level-crossing method can then be applied to 
measure the coupling between processes. 

Various techniques have been developed for studying 
the coupling between processes, e.g. the cross-correlation 
and cross-spectral density [18l |l9|, random matrix 
theory |20j- 24|] , the detrended cross-correlation analysis 
(DXA) 00, partial-DXA 01, coupled-DXA 0, 
and the detrending moving-average cross-correlation 
analysis [0. These methods miss the contribution of 
local effects due to the fact that the act of averaging 
constitutes their back bone, bringing up the idea of 
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using a method independent of averaging which only 
counts on local effects. This is where the quest for the 
generalized level crossing is intensified. In the present 
study we develop the level-crossing method to account 
for multiple stochastic processes. We find an analytic 
result for a multiple process whose components are 
independent Gaussian white noises. This result plays 
the role of a criterion of a fully uncorrelated multiple 
process. By comparing the analytic result with that of a 
coupled multiple discrete-time process, the efficiency of 
the coupling in a multiple process is obtained. 


B. Background 


Level crossing has been mainly implemented when 
dealing with continuous-time stationary stochastic 
processes, especially the famous class of Gaussian 
processes. In this subsection, we glance over the 
level-crossing method for a continuous-time stationary 
process A = {X t ,t > 0}. An upcrossing of the level 
x at time to is bound to occur if we have X t < x in 
(to — 5, to) and X t > x in (to,to + S), where S is the 
neighborhood radius. The number of upcrossings of the 
level x by X over the interval [0, t] is denoted by 7V t + (x). 
As A t + ( x) is itself a random variable, the mean number of 
upcrossings shown by ( N ( + (x)) is the quantity of interest. 
The general form of (7V t + (x)) for the continuous-time 
stationary process X is proved to be given by Q 

pOO 

(Nf(x))=t dz zp(x, z), (1) 

Jo 


where p(x,z) is the joint probability density of (X, X) 
containing all information about the level-crossing 
characteristics of a stationary process X. Note that the 
parameter t on the right-hand side of this equality comes 
from the stationarity of X. 

We turn our attention now to the case of stationary 
Gaussian processes. For a standardized stationary 
Gaussian process with correlation function r(r) = 
(XtXt+ T ) 3l|, the mean number of upcrossings of X is 
given by [31 


(N+(x)) = ^\ /2 exp y) , (2) 

where A 2 is the second moment of the spectral density 
function /(A) = r(r)e _lAr dr/2 t:. The Gaussian 
white noise is a special case in the family of stationary 
Gaussian processes whose correlation function is given by 
the Dirac delta function, <5(r). It should be noted that 
for a Gaussian white noise, Eq. (??) diverges because its 
second spectral moment A 2 is infinity. 

Numerical studies have shown that Eq. (??), which 
is for continuous-time stationary Gaussian processes, 
is inconsistent with the level-crossing results of 
discrete-time stationary Gaussian processes. Therefore, 


studying the level-crossing method for discrete-time 
processes could help avoid this discrepancy. The two 
following sections are devoted to this issue. In the fourth 
section, we create a criterion for measuring the coupling 
in multiple processes. Note that in this paper since we 
only deal with discrete-time processes, for simplicity the 
term “stochastic process” is used instead of the term 
“discrete-time stochastic process”, unless for emphasizing 
the type of process. 


II. LEVEL CROSSING FOR DISCRETE-TIME 
STOCHASTIC PROCESSES 


Consider a discrete-time stochastic process represented 
by {A} = X(t 2 ), • • •, X(t n )}. As stated earlier, 

the intention is to find the number of upcrossings, N+ (x), 
of a typical level, x, by the process {X} in the time 
period between t\ and t n . Mathematically speaking, an 
upcrossing of the level x at time tj is when we have 
X(ti- 1 ) < x and A(fj+i) > x. Since {X} is a stochastic 
process, 7V+(x) would be a random variable. Therefore, 
its ensemble average which is denoted by (7V+ (x)) is the 
quantity of interest. In case of a stationary stochastic 
process, (N+(x)) would become proportional to n with 
a proportionality constant represented by z/+. Note that 
is the average frequency of the upcrossings of the 
level x, where its inverse t x = l/v x gives the waiting 
time expressed as the average time expected for two 
consecutive upcrossings by {A}. The average frequency 
u+ for a stationary process {X} is given by 

/ X pOO 

dx 1 / dx 2 p x (x i,x 2 ), (3) 

-OO J X 


where px(x i,X 2 ) is the two-points joint probability 
density of the process {A} for two successive points Xi 
and X 2 . The right-hand side of Eq. (??) is nothing but 
the occurrence probability of an upcrossing of the level x 
by two successive points of {A}. 

For a standardized stationary Gaussian process, we 
have 


Px(xi,x 2 ) 


1 

27TA/1 - p 2 


exp 


x\ — 2px\x 2 + x| j 

2(1 -p 2 ) y 


(4) 

where p = (x\X 2 ) is the correlation between two 

successive points xi and X 2 . Substituting this equation 
into Eq. (??) gives 


vt = 


2 y^ J- 00 l 

(x/V2(l - p 2 )- pu 


1 — erf I 


(5) 


where 


erf(x) 



dt e 


( 6 ) 
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is the error function. It is important to note that the 
level-crossing result for discrete-time stationary Gaussian 
processes in Eq. (??) obviously differs from the result 
for continuous-time stationary Gaussian processes in 
Eq. (??). In the special case of no correlation, p = 0, 
i.e. when {X} is Gaussian white noise, Eq. (??) could 
be simply written as 


vt = 


1 

4 


1 - erf 


2 



(7) 


The average frequency in Eq. (??) for several values 
of the correlation p is depicted in Fig. [TJ The curve for 
p = 0 in Fig. |T| is for the case of Gaussian white noise. 
As seen in Fig. |T] the average number of upcrossings 
of the level x increases by decreasing the correlation p. 
This observation could be justified by the fact that a 
negatively (positively) correlated process, p < 0 (p > 0) 
fluctuates more (less) than a white noise (p = 0) which 
in turn leads to increasing (decreasing) the chance of 
upcrossing. It is noteworthy to state that the analytic 
description [Eq. (??)] provided in this section confirms 
the numerical results obtained for fractional Gaussian 
noises using the level-crossing method, see Ref. {32| . 


III. LEVEL CROSSING FOR DISCRETE-TIME 
VECTOR STOCHASTIC PROCESSES 


The level-crossing method proposed and implemented 
in studies prior to this work were basically able to 
analyze just a single continuous-time stochastic process 
or, at best, to analyze multiple processes in an individual 
manner (la . [33i ]. These are not options in this work, 
due to the fact that processes are not always individual 
or independent of each other. To comply with the 
aims of the present study where the coupling of the 
stochastic processes is not neglected, a vector stochastic 
process needs to be considered |Hj], and hence a new 
method needs to be provided. As such, we develop the 
level-crossing technique so we can call it the generalized 
level-crossing method. 

Mathematically speaking, a vector stochastic process 
is considered a d-dimensional path represented by {X} = 
{X(ti),X(t 2 ),... ,X(i n )}; where we have X(tj) = 
(Xi (ti), X 2 (ti),... ,Xd{ti)). The index d (> 2) indicates 
the number of stochastic processes being handled at the 
same time instance. To state clearer, when dealing with 
two simultaneous stochastic processes, e.g. the price 
return fluctuations of oil and gold, the parameter d 
would be equal to 2, which is due to the fact that two 
markets are being considered. We propose our apparatus 
for extracting statistical information on the path of a 
vector stochastic process {X}. This is what we call the 
generalized level crossing for a vector stochastic process 
in the d dimension. This method is founded on the 
combination of two concepts; radial and angular level 
crossings. 



Figure 1. The average frequency of the upcrossings of the 
level x for various values of correlation p in a discrete-time 
stationary Gaussian process. 


A. Radial level crossing 


Consider a d-dimensional vector stochastic process, 
{X}, and a d-dimensional sphere of radius R centered 
at the origin. In typical level-crossing methods, a level 
is crossed by a stochastic process. In the model studied 
here, the vector stochastic process ({X}) passes through 
the surface of a d-dimensional sphere, see Fig. UK a) for 
a three dimensional illustration. Now we define the 
radial level crossing for a vector stochastic process in 
a similar manner to the definition of the level crossing 
for a stochastic process. In the time period from t\ to 
t n the desired path {X} outcrosses our sphere N+(R) 
times. As in the one-dimensional case, an outcrossing of 
a d-dimensional sphere of radius R at time ti is when 
the conditions X(ti- 1 ) < R and X(ti+i) > R are 
complied. Note that X is the distance of the point X 
from the origin. Since {X} is a stochastic path, 7V+ ( R ) 
would posses a random behavior. Therefore its ensemble 
average, which is denoted by (N+(R)), would be our 
desired parameter. The average number of outcrossings 
depends on two facts: first, on the radial distribution of 
the points on the path {X} and, second, on the radial 
correlation of {X}. In the case of a stationary vector 
stochastic process, (N+(R)) would become proportional 
to n with a proportionality constant z/J [35l |. Since, 
z/^ is the average frequency of the surface (with radius 
R) outcrossings, its inverse represented by r# = 1/z/j^ 
gives the radial waiting time between two successive 
outcrossings. 

Quite similar to the one-dimensional case of the 
previous section, the average frequency v~^ in the d 
dimension for a stationary vector stochastic process {X} 
is given by 



where px(® 1 ,^ 2 ) is the two-point joint probability 
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density of the vector stochastic process {X} for the two 
successive points X\ and X 2 . So, the right-hand side of 
Eq. (??) is the occurrence probability for an outcrossing 
of the surface of the d-dimensional sphere with radius R 
by two successive points of {X}. Due to the importance 
of Gaussian processes we consider {X} to be a stationary 
vector stochastic process with standardized Gaussian 
distribution. In this case, we have 

iFvsr xp H pTj ^ 1 ' p }' (9) 

where P = (xi,x 2 ) t is a column vector containing 
coordinates of the two successive points, the superscript 
“T” indicates the transpose operation, and K = ( PP T ) is 
the covariance matrix. All information about correlations 
between the two points x\ and x 2 is placed within the 
covariance matrix K. In spite of having px.(xi,x 2 ) for 
Gaussian vector stochastic processes, the double integral 
of Eq. (??) cannot be further simplified unless it is for 
the special case of no correlation, where we refer to it in 
the next section. 

Although in this stage the radial level crossing 
is understood, valuable information from the vector 
stochastic process still could not be extracted, since the 
radial level crossing tells us nothing about the angular 
behavior of the vector stochastic process. This brings 
need for the presentation of the angular level crossing. 


B. Angular level crossing 

Consider one of the Cartesian axis Xi where, at the 
same time, is the axis of a cone, see Fig. d|b) for a 
three dimensional illustration. The cone apex is located 
at the origin of the coordinate system. The cone is 
specified by its apex angle, a, that takes values between 
0 to 7 r. Now the average number of outcrossings that 
path {X} experiences through the side surface of the 
cone represented by (N+(a)), is what we are going to 
count. Note that the average number of outcrossings 
depends on two facts. Firstly, on the angular distribution 
of the points on the path {X}, and secondly, on the 
angular correlation of {X}. In case of a stationary vector 
stochastic process, (_/V+(a)}, would become proportional 
to n with a proportionality constant Since, z/j is the 
average frequency of the outcrossings through the side 
surface of the cone, its inverse represented by r a = 1 /zz+ 
gives us the angular waiting time. The angular waiting 
time T a is the the average time between two consecutive 
side- surface outcrossings by {X}. 

As in the radial case, the average frequency 1 /+ for a 
stationary vector stochastic process {X} is given by 

= / d d x 1 f d d x 2 Px(x 1 ,x 2 ), (10) 

where px.(xi,x 2 ) is the two-point joint probability 
density of the vector stochastic process {X}. Angles 


(a) (b) 

x 2 x 2 




Figure 2. A schematic illustration of the radial (a) and 
angular (b) level crossings. Panel (a): A cut of a whole sphere 
limited to the positive axis of the Cartesian coordinates. The 
dashed-solid line sketches a path taken by a vector stochastic 
process named {X}. Panel (b): A cone where its axis overlaps 
with one of the Cartesian axes. Note that for both panels the 
black points are created when the path crosses the surface of 
the shapes, and the dashed and solid parts of the line indicate 
whether its inside or outside the shapes. 

and fl 2 in Eq. (??) are made by the vectors X\ and 
x 2 with the Xi axis Respectively. When {X} represents 
a stationary Gaussian vector stochastic process, the 
same explanations expressed for the radial level crossing 
applies. 


IV. THE CREATION OF A CRITERION 

However, we are still not there yet. The reason is 
that when implementing Eqs. (??) and (??) for radial 
and angular level crossings, the results are not conclusive 
by themselves. This is due to the fact that the results 
should first be valued. In other words they need to be 
compared with some sort of a criterion to provide a basis 
for the most suitable and applicable conclusions. Usually 
the best criterion that would work as a measure must 
not be biased. Therefore, in the context of the present 
study, the criterion is selected to be an uncorrelated 
Gaussian process. The intention is to obtain an analytic 
solution for the frequencies and z/+ of a vector 
stochastic process in d dimensions whose components are 
independent Gaussian white noises. 

To obtain the criterion, consider a vector stochastic 
process consisting of d independent Gaussian white noises 
represented by {X„}. The joint probability distribution 
for two successive points x\ and x 2 of the path {X mn } 
is given by Eq. (??) in which the covariance matrix I\ 
is equal to the identity matrix I. The joint probability 
distribution px(®i 7 * 2 ) is then reduced to 

Px{xi,x 2 ) =p(x 1 )p(x 2 ), (11) 


p(x) = — 7 = —exp ( — -x.x T ). ( 12 ) 

K {V2^) d V 2 J 


with 
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Equation (??) means that there is no correlation between 
the two points leaving them independent of each other. 
Using this equation, the double integral of Eq. (??) is 
split into the product of two univariate integrals. So, for 
the vector stochastic process {X„, n }, the radial frequency 
is given by 

U+ = P in P ouU (13) 

where 


origin, the axis overlaps with the X* axis, and the apex 
angle is denoted by a. By substituting the probability 
density function of Eq. (??) into Eq. (??) and integrating 
over the solid angle of the cone, the probability Pi n is 
obtained as follows: 

{ a/27r d = 2 

( 20 ) 

sin 2 (a/4) d = 3, 


Pin = / d d x 1 p{x i), (14) 

J xi <R 

and 

Pout = [ d d x 2 p{x 2 ). (15) 

J X2 >R 

As seen from the domain of integration in Eqs. (??) and 
(??), Pi n represents the probability of the first point to 
be inside the sphere of radius R centered at the origin, 
and P ou t is the probability of the second point to be 
outside the very sphere. Since every point on the path 
is either inside or outside the sphere, the summation of 
the probabilities P. m and P ou t is equal to unity. Using 
Eq. (??) the probability ( Pi n ) is obtained as follows: 

o(l-d/2) rR 

Pin = (d/2 — 1)1 J 0 tex-'T'l . (16) 

The integration of Eq. (??) can be calculated using 
the method of integration by parts. This enables an 
analytical solution for the positive frequencies, where by 
substituting the expressions for Pi n and P ou t = 1 — Pin 
from Eq. (??) into Eq. (??), is obtained. 

Due to the isotropic characteristic of the path {X„,„}, 
there is no preferred direction. Therefore, we could 
obtain i/+ in Eq. (??) for an arbitrary axis, Xi , and say 
that this result would be the same for all other directions. 
The angle between a typical point x on the curve {X^} 
and the x^ axis is denoted by d. Since the two successive 
points Xi and x 2 on the path {X m „} are independent 
random vectors, the angles d\ and d 2 are independent 
random variables. This independence enables us to write 
the similar equation as Eq. (??) for the case of So 
the angular frequency is given by 

V+ = P in P ouU ( 17 ) 

with 

Pin = [ d d x 1 p(x i), (18) 

J i?i <a/2 

and 

Pout = [ d d x 2 p{x 2 ). (19) 

j $2 >ct /2 

Note that, in this case, P ln ( P ou t ) is the probability that a 
point resides inside (outside) a cone with the apex on the 


for the lower dimensions d = 2,3, and 


p 

rr 



r a/2 


dd (sin$) 


d—2 


( 21 ) 

for the higher dimensions d > 3, where T(x) is the 
gamma function. This enable us to have an analytical 
solution for the positive frequency, = Pi n P ou t , which 
is obtained by the expressions for P in and the equality 
Pout = 1 — Pin- Note that for extracting Eq. (??) we 
used the hyperspherical coordinate system which is the 
generalization of the spherical coordinate system to the 
dimensions higher than 3. 


The interest here is to browse the collective behavior 
of a multiple process by crossing a specific level. In other 
words, we show how coupling is featured in the context 
of the present study. Toward this end, we investigate 
the generalized level-crossing method in two dimensions. 
The positive frequencies for the criteria in two dimensions 
is obtained by substituting the expressions for and 
Pout = 1 — Pin from Eqs. (??) and (??) into Eqs. (??) 
and (??), 


v+(R)=e~ R2 ' 2 ( K l-e- fl2 / 2 ), (22) 

and 



where the zero index is used to emphasize that these 
frequencies are for the uncoupled case and distinguish 
them from the coupled case. Now, consider a 
vector stochastic process {X} = {(xi, y±), (x 2 , y 2 ),... } 
consisting of two standardized Gaussian white noises 
{xi, x 2 ,... } and {j/i, y 2 ,... } with a Gaussian coupling 
as 


C{n) = {Xiy l+n ) = Aexp , (24) 

where (.) denotes the ensemble average, U A V represents 
the amplitude of coupling, and is the correlation 
length. In this case, there exists no analytical expression 
for the average frequencies and z/+. Therefore they 
could only be computed by the numerical integrations of 
Eqs.(??) and (??). Note that Eq. (??) determines the 
elements of the covariance matrix, K : in Eq. (??). 
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In order to show the deviation between the frequencies 
{ v Ri v a) °f coupled white noises and (uq(R), Vq (a)) of 
uncoupled white noises, we introduce the following ratios 


= Vr ~ Vp (R) 
v+(R) 


(25) 


and 


D 


ang — 


~ "0 («) 

Vo( a ) 


(26) 


Figure [3] shows these deviation ratios for radial [Fig. [31(a)] 
and angular [Fig. EKb)] crossings for the Gaussian 
coupling in Eq. (??) with the amplitude A = 1/3 
and correlation lengths £ = 1,2,3,4. Notice that 
the deviation is directly proportional to the correlation 
length, where a higher correlation length gives a curve 
with a bigger deviation. This is true for both radial and 
angular level crossings. Another important conclusion 
made from Figs. Ha) and[3|b) comes from the size of the 
correlation length £. The fact of the matter is that this 
itself is a standing point for the promise of the present 
study. In other words, although the correlation length is 
very small, there exists a pronounced deviation between 
the coupled and uncoupled cases. 


V. CONCLUSIONS 


this statement. The generalized level-crossing method 
consists of two working concepts, namely radial and 
angular level crossing. These new concepts enable us 
to study the coupling between processes. In order 
to pronounce the efficiency of the coupling between 
the components of a multiple process, we introduced 
the state-of-the-art criterion which is the benchmark of 
multiple processes with no coupling. We derived analytic 
results for the radial and angular level-crossing regarding 
this criterion, see Eqs. (??) and (??). In order to 
evaluate the criterion, the generalized level-crossing has 
been studied in two dimensions. The results show the 
sensitivity of the radial and angular solutions to slight 
couplings with short correlation lengths. 


(a) 



All processes in nature, although they may seem 
to be continuous, are actually not. As a matter of 
fact, all processes are discrete, time wise. In other 
words, a seemingly continuous process is a look from 
far at that process. Now by looking closer and closer, 
its discreteness becomes observable. Toward that end, 
we extended the level-crossing method to the realm of 
discrete-time stochastic processes, which enabled filling 
two of the important gaps in this method. First, we 
obtained an analytical expression for the level crossing 
of a discrete-time stationary Gaussian process, see 
Eq. (??). The generality of the expression is due to 
the fact that it is derived for a Gaussian process with 
an arbitrary correlation. Second, we developed the 
level-crossing method to enable simultaneous analysis 
of several discrete-time processes. The reason for 
going this way is due to the fact that processes are 
not exactly independent. Our analytic modeling and 
hence solutions contribute towards better understanding 


(b) 



Figure 3. Sensitivity of the two dimensional level-crossing 
analysis on coupling with short correlation lengths. Panels (a) 
and (b) are the result of comparison between two uncoupled 
white noises and two coupled white noises, respectively, in the 
radial and angular level-crossing regimes. 
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